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Abstract. We present a possible interpretation for the increase of the angular diameter of the Mira variables o Get, R Leo, 
and X Gyg from the K band to the 1 1 yum region revealed by the recent interferometric observations using narrow bandpasses 
where no salient spectral feature is present (Weiner et al. 2003a 2003b i. A simple two-layer model consisting of hot and cool 
H2O layers for the warm water vapor envelope, whose presence in Mira variables was revealed by previous spectroscopic 
observations, can reproduce the angular diameters observed with Infrared Spatial Interferometer as well as the high-resolution 
TEXES spectra obtained in the 1 1 yum region. The warm water vapor layers are optically thick in the lines, and therefore, strong 
absorption due to H2O can be expected from such a dense water vapor envelope. However, the absorption lines are filled in by 
emission from the extended part of the envelope, and this results in the high-resolution 1 1 yum spectra which exhibit only weak, 
fine spectral features, masking the spectroscopic evidences of the dense, warm water vapor envelope. On the other hand, the 
presence of the warm water vapor envelope manifests itself as the larger angular diameters in the 1 1 jim region as compared 
to those measured in the near-infrared. Furthermore, comparison of the visibilities predicted in the near-infrared (1.5-3.8 yum) 
with observational results available in the literature demonstrates that our two-layer model for the warm water vapor envelope 
can also reproduce the observed near-infrared visibilities and angular diameters, and suggests that the wavelength dependence 
of the angular size of Mira variables in the infrared largely reflects the H2O opacity. The radii of the hot HtG layers in the three 
Mira variables are derived to be 1.5-1.7 Rt, with temperatures of 1800-2000 K and H2O column densities of (1-5) x 10^' cm"^, 
while the radii of the cool H2G layers are derived to be 2.2-2.5 R^, with temperatures of 1200-1400 K and H2O column densities 
of (1-7) X 10^' cm-2. 

Key words, infrared: stars - molecular processes - techniques: interferometric - stars: late-type - stars: AGB and post-AGB - 
stars: individual: o Get, R heo,x Cyg 



K*" 1 . Introduction Although these Mira variables exhibit dust emission from 



the circumstellar envelope, the observed increase of the angu- 
An increase of the angular size from the near-infrared toward j^ diameter from the near-infrared to the mid-infrared cannot 



d longer wavelengths appears to be a common phenomenon ob- ^^^^^^ ^^ attributed to the dust shell, as discussed in detail by 
served m various classes of cool luminous stars s uch as M wHT03a. The radii of the inner boundary of the circumstellar 
(super)giants and Mira variables. Mennesson et al. (l002 1 re- ^^^j ^jj^jj ^^ estimated to be 2-4 R^ for o Cet and R Leo, and 
vealed that the umform disk (UD) diameters of semiregu- ^20 R. for ;r Cyg (Danchi et al.|IMI Schuller et aLEUPl. 
lar M giants as well as Mira variables measured in the L' ^^^^ ^^^j^ extended dust shells are completely resolved with the 
band are 20-100% larger than t hose me asured in the K' band, baselines of 20-56 m used in the ISI observations by WOO, 
Recently, Weiner et al. (2000 !2003al and 2003b. hereafter wHT03a, and WHT03b. In other words, the visibility compo- 
WOO, WHT03a, and WHT03b, respectively) have measured the ^^^^ resulting from the dust shell is nearly zero at these baseUne 
angular size at 11 ^/m for the M supergiants a On and a Her as j^^g^j^^^ ^„j tj^^ g^p^^j ^f ^^^ ^^^t ^jj^H -^ ^^ j^^^^ ^j^^ ^^^^j ^■^_ 
well as the Mira variables o Cet, R Leo, and x Cyg, using the Ability by an amount equal to the fraction of flux originating in 
Infrared Spatial Interferometer (ISI) with a spectral bandwidth ^^e dust shell in the field of view. This eff^ect is already taken 
of 0. 17 cm-i . They found that the uniform disk diameters of -^^^ ^^^^^^^ -^ ^j^^ determination of die uniform disk diameters 
the two M supergiants measured at 1 1 ^m are -30% larger than ^y WOO, WHT03a, and WHT03b, and therefore, it is unlikely 
those measured in the K band, while the 1 1 fim uniform disk ^^at the extended dust shell causes the 1 1 /zm angular size to 
diameters of the three Mira variables are roughly twice as large ^pp^^^ ^^^^j. ^^^ j„ ^^e near-infrared. 
as those measured in the K band. 
Instead, the increase of the angular size from the K band to 
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et al. (I2002> conclude that the observed increase of the an- 
gular diameter from the K' band to the L' band can be ex- 
plained by an optically thin layer (r ^ 0.5) extending to ~3 /?* 
with a (pseudo)continuous opacity. Perrin et al. t2004j apply 
a similar model to interferometric observations of a Ori and 
a Her obtained in the K and L bands as well as in the 1 1 /im 
region, and conclude that the extended gaseous layer is opti- 
cally thin in the K and L bands {tk - 0.06 and ti - 0.026), 
but optically thick in the 11 jum region (rn^m = 2.33). Both 
Mennesson et al. ( 2002 1 and Perrin et al. l2004j suspect opac- 
ities due to molecular species such as H2O, CO, and SiO as 
the source of the (pseudo)continuous opacity. The presence of 
such a warm molecular envelope is akeady confirmed by anal- 
yses of infrared molecular spectra. Tsuji P978 1 interpreted in- 
frared spectral data of the M supergiant yu Cep and the Mira 
variable R Cas in terms of a warm water vapor envelope, 
and Tsuji ( I1988> further found evidence of the presence of 
the warm molecular envelope in the high-resolution spectra 
of M (super)giants. More recently, infrared spectra obtained 
with the Short Wavelength Spectrometer (SWS) onboard the 
Infrared Space Observatory (ISO) have provided ample infor- 
mation on physical parameters of the warm molecular envelope 
(e.g., Tsuji et al. [T9^ Yamamura et al. [T999l Cami et al. l20nol 
Matsuuraetal. l2002> . 

However, the high-resolution 1 1 /im spectra of the M su- 
pergiants a Ori and a Her as well as the Mira variables o Cet, 
R Leo, and x Cyg presented by WHT03a and WHT03b turn 
out to pose difficulty in interpreting the increase of the angu- 
lar diameter in terms of the extended, warm water vapor en- 
velope. These 1 1 /im spectra were obtained using the TEXES 
instrument mounted on the Infrared Telescope Facility tele- 
scope with a spectral resolution of ~10^ (Lacy et al. 2002 1. 
The high-resolution TEXES spectra of the above M supergiants 
and Mira variables reveal that no salient spectral features are 
present within the bandpasses used in the ISI observations, and 
therefore, these spectroscopic observations appear to be in con- 
tradiction with the interpretation that the dense, warm water va- 
por envelope is responsible for the increase of the angular size 
from the near-infrared to the 1 1 //m region. 

Ohnaka (.2004. hereafter Paper I) demonstrates that the 
warm water vapor envelope extending to L4-L5 R* with tem- 
peratures of ~2000 K and H2O column densities of the order 



of 10' cm can reproduce the observed increase of the an- 
gular diameter from the near-infrared to the 1 1 //m region in 
a Ori and a Her and, simultaneously, their featureless high- 
resolution spectra at 11 fiva. While significant absorption is ex- 
pected from such a warm water vapor envelope, the absorption 
lines are filled in by the emission of H2O lines from the ex- 
tended part of the envelope, which leads to spectra without any 
conspicuous spectral features as observed with the TEXES in- 
strument. On the other hand, the emission from the extended 
water vapor envelope causes the apparent size of the objects 
to be significantly larger than the photospheric size, in agree- 
ment with the interferometric observations with ISI presented 
by WHT03b. 

If the warm water vapor envelope is present even in early 
M supergiants such as a Ori and fi Cep and causes the stel- 
lar angular size in the mid-infrared to appear larger than that 
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Fig. 1. The model used in the present work. The photosphere 
is approximated by a blackbody of Teft with a radius /?*. The 
column densities in the radial direction (A^hot and A^cooi) and the 
temperatures of H2O gas (Thot and Tcooi) as well as the radii of 
the water vapor layers (7?hot and /?cooi) are the input parameters. 

measured in the near-infrared, it can also be the case for Mira 
variables, for which water vapor is known to be present in the 
photosphere and we can expect even denser water vapor en- 
velopes than in early M supergiants. In the present paper, we at- 
tempt to interpret the increase of the angular diameter from the 
near-infrared to the 1 1 /im region as well as the high-resolution 
1 1 //m TEXES spectra observed for the Mira variables o Cet, 
R Leo, and x Cyg in terms of the warm water vapor envelope. 

2. Modeling of the warm water vapor envelope 

In our model, the central star is surrounded by a warm water va- 
por envelope, which consists of two layers with different tem- 
peratures Thot and Tcooh as depicted in Fig.^ The central star is 
approximated by a blackbody of effective temperature T^s with 
a radius /?*. The temperature and the density are assumed to be 
constant in these hot and cool layers, which extend to /?hot and 
Rcoo\, respectively. We denote the column density of H2O in the 
radial direction in the hot and cool layers as A^hot and Ncooh re- 
spectively. As very simple it may seem, such a two-layer model 
turned out to be successful in interpreting infrared molecular 
spectra obtained with ISO/SWS. For example, Yamamura et al. 
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<1999i and Matsuura et al. ('2002^ could reproduce H2O spec- 
tra of Mira variables obtained with ISO/S WS using a two-layer 
model, although their two-layer model consists of a stack of 
plane-parallel slabs, instead of spherical layers as adopted in 
the present work. 

We adopt an effective temperature of 3000 K for the Mira 
variables studied in the present work. While the effective tem- 
peratures of Mira variables are known to show temporal varia- 
tions with phase (e.g., van Belle et al. 1996 Perrin et al. 1999 
Woodruff et al. I2004> . the effect of the uncertainty of the ef- 
fective temperature adopted in our model turns out to be minor 
as we show below, because the warm water vapor envelope is 
quite optically thick in the wavelength regions discussed in the 
present work. 

We first calculate the line opacity due to H2O assuming a 
Gaussian profile. We adopt a velocity of 5 km s"' for the sum 
of the thermal velocity and the micro-turbulent velocity. The 
gf -values and the lower excitation potentials of H2O lines are 
calculated using the HITEMP line list (Rothmanjl997 1 and the 



list compiled by Partridge & Schwenke J 19971 hereafter PS97). 
The HITEMP line list is an extension of the HITRAN database 
toward high temperatures expected in stellar atmospheres, and 
the HITEMP H2O line list used here is calculated for a tem- 
perature of 1000 K. As we demonstrate below, the temperature 
of the warm water vapor envelope can be as high as ~2000 K, 
and therefore, the HITEMP H2O line list for 1000 K may not 
be adequate for such a high temperature. The PS97 H2O line 
list is generated for a temperature of 4000 K and appears to be 
the most extensive line list of H2O up to date. We use both line 
lists for our calculation to study possible effects of the differ- 
ence between the two H2O line lists. 

The energy level populations of H2O are calculated assum- 
ing local thermodynamic al equilibrium (LTE). We examine the 
validity of LTE in the temperature and density ranges relevant 
for the present work, using the order-of-magnitude estimates 
of collisional and radiative de-excitation rates, as adopted in 
Paper I. The collisional de-excitation rate Cui is given by Cui ~ 
A^cTuiVjei, where A^ is the density of the primary collision part- 
ner, (Tui is the cross section, which we approximate with the 
geometrical cross section, and Viei is the relative velocity be- 
tween the colUsion partner species and H2O molecules. As we 
will show below, for the hot layer of the warm water vapor 
envelope in o Cet, R Leo, and x Cyg, the column density of 
H2O is found to be of the order of 10^' cm"^ with tempera- 
tures of 1800-2000 K, and the geometrical thickness is derived 
to be ~0.5 Ri,, which translates into 1 x lO'"* cm with a stel- 
lar radius of 300 Rq assumed. The number density of H2O 
is then estimated to be 1 x 10^ cm""^. For temperatures and 
densities found in the hot layer of the warm water vapor en- 
velope, H atoms appear to be the primary collision partner of 
H2O molecules, and the ratio of the number density of H atoms 
to that of H2O molecules expected in chemical equilibrium is 
approximately 10"* for the relevant temperatures and densities. 
Therefore, the number density of H atoms in the hot layer is 
estimated to be ~ 1 x 10'^ cm"-'. With a geometrical cross sec- 
tion cTyi of 10"'^ cm^ and a relative velocity Vrei of 5 km s"' 
assumed, this number density of H atoms leads to a collisional 
de-excitation rate of ~ 500 s '. For the cool component of the 



water vapor envelope, we derive temperatures of 1200-1400 K, 
geometrical thicknesses of ~ 1 .5 /?*, and H2O column densities 
of the order of 10^' cm"^, which leads to a number density of 
H2O of ~ 3 X 10^ cm"^. At these temperatures and densities, H2 
molecules are the main collision partner of H2O molecules, and 
the ratio of the number density of H2 molecules to that of H2O 
molecules expected in chemical equilibrium is ~ 3 x 10^, which 
leads to a collisional de-excitation rate of ~ 50 s ' . On the other 
hand, the rate of spontaneous emission can be estimated from 
the Einstein coefficients A^x. For the H2O molecule, the ranges 
of Aui ai-e approximately ~10"'^-3 x 10^ s"', ~10"^-3 x 10^ s"\ 
and % 10 s' in the 1 1 yum region, K band, and L' band, respec- 
tively. Therefore, for these wavelength regions we will discuss 
below, the assumption of LTE is valid for weak and moder- 
ately strong H2O lines both in the hot and cool layers of the 
warm water vapor envelope, while non-LTE effects may not be 
negligible for very strong lines. However, a comprehensive cal- 
culation with non-LTE effects taken into account is beyond the 
scope of the present paper, and we assume LTE for the H2O 
lines considered here. 

Once the line opacity has been obtained, the monochro- 
matic intensity profile as well as the spectrum (emergent flux) 
can be calculated at an appropriate wavenumber interval, as 
described in Paper I. In order to compare with the observed 
TEXES spectra of the three Mira variables, synthetic spectra 
are spectrally convolved with a Gaussian profile which repre- 
sents the effects of the spectral resolution of the instrument as 
well as of the macro-turbulent velocity in the atmosphere of 
the Mira variables. The spectral resolution of the TEXES in- 
strument is 10'', which translates into an instrumental broad- 
ening of 3 km s '. The macro-turbulent velocity in the atmo- 
sphere of Mira variables is not well known. For non-Mira M 
giants, however, analyses of high-resolution spectra are avail- 
able, and they may be used as a guide to estimate the macro- 
turbulent velocity in the atmosphere of Mira variables. Tsuji 
( 1986 1 analyzed the CO first-overtone bands of non-Mira K and 
M giants and obtained macro-turbulent velocities ranging from 
~1 kms"' to ~4 kms '. In the present work, we tentatively 
assume a macro-turbulent velocity of 3 km s"' for the three 
Miras, and the FWHM of the Gaussian profile to represent the 
spectral resolution of the instrument and the macro-turbulent 
velocity is calculated as V3^ + 3^ - 4.24 kms"'. The spec- 
trum convolved with this Gaussian profile and normalized, F^, 
is then diluted as follows, taking the effect of continuous dust 
emission from the circumstellar dust shefl into account: 

^w " "^ = (1 - fdust)Fa, + /dust , (1) 

where /^^iiuted j^ ^^^ ^^^^ spectrum and /dust is the fraction of 
the flux contribution of the circumstellar dust shell. 

We calculate the monochromatic visibility from the 
monochromatic intensity profile at each wavenumber using 
the Hankel transform, and then the monochromatic visibility 
is spectrally convolved with an appropriate response function 
which represents the spectral resolution of the interferometric 
observation at issue. For the ISI observations we will discuss 
below, WHT03a and WHT03b used narrow bandpasses with a 
width of 0.17 cm"'. Therefore, the spectral response function 
is assumed to be a top-hat function with a width of 0. 17 cm"'. 



K. Ohnaka: Warm water vapor envelope in Mira variables 



We compare synthetic spectra and uniform disk diameters 
with the TEXES spectra as well as the 1 1 jjm uniform disk di- 
ameters obtained by WHT03a and WHT03b, searching for a 
parameter set which can best reproduce these observational re- 
sults. The ranges of the input parameters are as follows: Thot 
(K) = 1600 ... 2000 with AThot = 100 K, Teooi (K) = 1000 ... 
1400 with Areooi = 100 K, T^hot (R*) = 1-3, 1.5, 1.7, and 2.0, 
and /?cooi (R*) - 2.0, 2.2, 2.5, 2.7, and 3.0 with a condition 
^hot < ^cooi- The grid for A^hot (cm^^ ) and A^cooi (cm^^ ) is 10'^, 
lO^", 5 X lO^o, 10^1, 3 X lO^i, 5 X 10^', 7 x lO^', and 10^1 
The uncertainties of the derived temperatures and radii of the 
H2O layers are estimated to be ±100 K and ±0.2 7?*, respec- 
tively, while the uncertainties of the H2O column densities are 
approximately a factor of 2. The adoption of an effective tem- 
perature of 2800 K, instead of 3000 K, leads to no noticeable 
changes in the spectra and visibilities discussed in the present 
work. In the next section, we discuss comparison between ob- 
servation and model for each object. 



3. Comparison with tlie observed spectra and 
angular diameters 

3.1.oCet 

3.1.1. 11 //m spectrum 

For o Cet, we find that the 1 1 //m spectrum observed at phase 
0.36, which is presented in WHT03a, as well as the 1 1 //m uni- 
form disk diameters obtained by WOO and WHT03a can be best 
reproduced by a model with Thot = 1800 K, /?hot = 1-5 R*, 
A^hot = 5 X lO^i cm-2, r.ooi = 1400 K, /Jcooi = 2.2 R^, and 
A^cooi = 1 X 10^' cm"~. Figure 121 shows the synthetic spectra in 
the 1 1 /im region predicted by the best-fit model for o Cet, to- 
gether with the TEXES spectrum obtained at phase 0.36, which 
was read off Fig. 2 in WHT03a. The synthetic spectrum cal- 
culated with the HITEMP H2O line list is represented by the 
green solid line, while that calculated with the PS97 H2O line 
list is represented by the blue solid line. The synthetic spec- 
tra are redshifted by 83 km s"' to match the observation (see 
WHT03a). WHT03a derived the fraction of the flux contribu- 
tion of the stellar disk of o Cet ranging from 0.24 to 0.47 at 
1 1 //m, which corresponds to the fraction of the flux contri- 
bution of the dust shell between 0.76 and 0.53. We adopted 
/dust = 0.7 in our calculation of the 1 1 jjm region for o Cet. 
While the synthetic spectrum calculated with the PS97 line list 
shows H2O features which are not seen in the spectrum calcu- 
lated with the HITEMP database, the discrepancy between the 
two synthetic spectra is not very serious or systematic for the 
temperatures and densities relevant for the present calculation, 
given the uncertainties of the Une positions of both line lists. 
For example, Ryde et al. ("2002"! found that the uncertainty of 
the positions of H2O pure rotation lines around 1 2 //m predicted 
by PS97 is ~0.05 cm '. On the other hand, the uncertainty of 
H2O line positions in the HITEMP database is between ~0.1 
and ~1.0 cm '. Figure |2] shows that the two synthetic spec- 
tra can reproduce the approximate depths and heights of the 
observed fine spectral features and, in particular, the observed 
spectra free from any conspicuous spectral features in the band- 





III 


II Hill 


1 III nil 1 1 nil 1 ri 


II II Id III iiiiiiii II iKiiiiiiii 1 


I'l iiiiiir 


II II mill 


iiiiriiiii 1 iiriiiiiii II mil i iiii 


1 IIII mil 


iiiir iiiiiiiiiii 


1.2 








f. 










- 


"^3 ■> 


^ 


r^ 


-Y^ 


^S?^W(30^ 


A 


^ 


»J'=?^^Q%J^ 


>Q! 


25??2>« 


o.s 


l\ 1 


II II 


II 1 1 


II III 1 1 mil II Mil 


iii'-iiii 1 
1, 


'ill 1 

1 

1 , 


1 llllll III III 1 1 II lllllllll 


1 III II IIII 


III 1 mil- 



918 

II liiiiiiiii mil inn I'liiiiiiiii II 



918.5 

0) (om"^) 

II III II 1 1 ll'l II I IIII II IIIIIIII llllll I III IIII 



919 

I III III 




902 
<. (cm-') 



mil 1 1 nil 


1 III I'l II II II III' 


III 1 Kill IIII IIIKIIIIIIII 


III' null 


iniiiii II 


lllllllll 


1 r 1 iiiiriiiiii 


ii'iiiiiiiiiiiiriiiiiiiiiiiir II II 
i 


-jS"^^ 


Y^^ 


7*0;^/*^^^ 


^«^s<L 


v^ 


1 

f 


1-Ow 




^"k^^ 


r^i 


HI 1 II 1 1 II 1 


II III III II 1 


1 II III II 1 nil 


IIII llllll llllll '1 

1 

, 1 


1 II 1 1 


l1 


1 Hill Mill III 


1 1 II 1 1 iiiiiiiiii II III ii-f 









0, 


(cm-') 














mil riiiiiiim I'm i iiiiii'i 


1 1 in'iii mil II lllllllll 


111 III nil 


III llllll 


IIII ll 


mil II (III II 


iiiimiiiim I'l III 


iiiiiiiiiiiii 


II III 1 


1.2 


-\J6^ 


^^x^^ 


y^^ 


4^ 


^ 




. 


>^ 


r< 


^- 


^ 


Q(^\ 


0,8 


hi 1 II 1 mil 1 1 iiiiiii 


1 II III III 1 II 





1 1 III' 1 1 
1 
1 . 


IIII 


1 1 II II 


II 1 II II II 


1 IIIIIIIIII 


1 in 



895 895.5 

„ (om->) 



Fig. 2. Spectra of o Cet in the 1 1 yum region at phase 0.36. The 
red dots represent the observed spectrum of o Cet presented in 
WHT03a, while the green and blue solid lines represent the 
calculated spectra using the HITEMP database and the PS97 
line list, respectively. The parameters of the best-fit model for 

Cet are That = 1800 K, R^ot = 1.5 R^, A^hot = 5 x lO^' cm^^, 
r.ooi = 1400 K, T^cooi = 2.2 R^, and A^^ooi = 1 x 10^' cm^l The 
synthetic spectra are convolved with a Gaussian with a FWHM 
of 0.013 cm ' to account for the effects of the instrument as 
well as of the macro-turbulent velocity, and are redshifted by 
83 km s ' to match the observation. The dashed lines repre- 
sent the bandpasses used in the ISI observations. The positions 
of the H2O lines whose intensity at 2000 K is stronger than 

1 X 10"^^ cm molecule"' are marked with upper ticks (PS97 
line list) and lower ticks (HITEMP database). These line posi- 
tions are also redshifted by 83 km s ' with respect to the rest 
wavenumber. 



passes at around ~896.83 cm"' and ~895. 1 cm"' (dashed lines) 
used for the ISI observations. In these bandpasses, a number 
of H2O lines are present, as marked by the ticks in Fig. |2l 
The continuum-like appearance of the spectra is a result of the 
filling-in of the absorption lines by emission from the extended 
part of the warm water vapor envelope, as we will discuss be- 
low. On the other hand, the match for the individual line profiles 
is not very good, which can partially be attributed to the uncer- 
tainties of the line positions predicted by both line lists, but it is 
also due to the simplicity of our model, which cannot account 
for complicated physical and chemical processes taking place 
in the outer atmosphere of Mira variables. For example, the 
strong H2O fines at 897.40, 897.45, 902.04, and 918.57 cm' 
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Fig. 3. a: 11 fim spectra predicted from the best-fit model for 
o Cet at the stellar disk center (//=!, top spectrum), near the 
limb of the warm H2O envelope (jj = 0.2, bottom spectrum), 
and between the disk center and the limb (// = 0.85, mid- 
dle spectrum). The continuum level is shown with the dashed 
line, and the HITEMP database is used in the calculation. The 
flux contribution of dust emission is not included in the spec- 
tra shown, b; The dotted line represents the spectrum (emer- 
gent flux) expected for the whole object (stellar disk + warm 
water vapor envelope) without the dilution effect due to dust 
emission. The solid line represents the final spectrum with a 
flux contribution of dust emission of 70% taken into account. 
c: Optical depth in the radial direction due to H2O lines in the 
same spectral region. The positions of the H2O lines are marked 
with ticks. The spectra and the flux shown in a and b as well as 
the optical depth shown in c are redshifted by 83 km s"' with 
respect to the rest wavenumber to match the observed spectrum 
shown in Fig.|2l but not convolved with the Gaussian represent- 
ing the instrumental eff'ect and the macro-turbulent velocity. 

show inverse P-Cygni profiles, which suggest an inflow with 
a velocity of 11 km s ' . Such a dynamical process is not taken 
into account in our model, and therefore, these line profiles can- 
not be well reproduced by our two-layer model. 

The continuum-like appearance of the spectrum of o Cet in 
the 1 1 /vm region can be explained by the filling-in effect due 
to H2O line emission from the extended part of the warm water 
vapor envelope. This effect is illustrated in Fig. |3 where the 
spectra expected at the disk center (// = 1, where fi = cosO 
with defined in Fig.Q, at a position between the disk center 
and the limb (// = 0.85), and near the limb (ji = 0.2) are shown 
in the top panel. The bottom panel of Fig. |3] shows the opti- 
cal depth of the water vapor envelope, and the figure reveals 



that the water vapor envelope is optically thick (r^ » 1-60). 
These spectra as well as the optical depth are calculated with 
the HITEMP database. The H2O lines appear as strong absorp- 
tion in the spectrum at the disk center, as the top spectrum in 
Fig. |3t shows. On the other hand, the spectrum near the limb 
(bottom spectrum in Fig.|^) shows emission at the positions of 
the H2O lines, and these emission lines originate in the cooler 
H2O layer The spectrum for a line of sight between the disk 
center and the limb (middle spectrum in Fig.|3^) is dominated 
by the H2O emission lines originating in the hotter layer over- 
laid by the H2O absorption lines originating in the cooler layer. 
Since the spectrum (emergent flux) of the object is obtained by 
integrating the intensity over the area that the object projects 
onto the plane of the sky, the absorption due to H2O lines is sig- 
nificantly filled in by emission from the limb, resulting in the 
spectrum showing only moderately strong emission features, 
as plotted with the dotted line in Fig.O). These emission fea- 
tures are further weakened by dilution due to continuous dust 
emission, as the final spectrum (solid fine) in Fig.|3j5 shows. 

The parameters derived above are similar to those derived 
by Yamamura et al. < 19991 1 for o Cet at phase 0.99: Thot - 
2000 K, 7?hot = 2.0 R^, Nhot = 3.0 x lO^' cm^^, Tcooi = 1400 K, 
Rcooi - 2.3 Ri,, and A^cooi - 3.0 x 10^" cm"^. However, it is 
not straightforward to discuss the phase dependence of these 
physical properties, not only because of the uncertainties of the 
derived parameters, but also because of the phase dependence 
of the stellar continuum radius. The analysis by Woodruff 
et al. tZOQ A") demonstrates that the stellar radius of o Cet in 
the 1 .04 jum continuum predicted by dynamical model atmo- 
spheres shows a sinusoidal variation from phase 0.0 to 0.5 with 
a maximum at around phase 0.25, while the relatively large er- 
rors of the observationally derived continuum radii resulting 
from the uncertainty of the distance prevented the authors from 
confirming this theoretical prediction. Interferometric observa- 
tions with spectral resolution high enough to resolve molecular 
bands would be useful for obtaining reliable stellar continuum 
diameters and studying temporal variations of physical proper- 
ties of the warm water vapor envelope. 



3.1 .2. 1 1 /im angular diameter 

We now calculate the angular diameter of o Cet predicted from 
the above model for the warm water vapor envelope. In order to 
compare with the observed uniform disk diameters, it is neces- 
sary to estimate the angular size of the star, which corresponds 



to Ri,. Woodruff et al. J2004t analyzed the A'-band interfero- 
metric data of o Cet obtained with VLTI/VINCI at phases from 
0.13 to 0.40. They derived the Rosseland radius as well as the 
continuum radius at 1.04 fim at each phase by fitting the ob- 
served visibilities with theoretical ones predicted by dynami- 
cal model atmospheres, and found that the continuum diame- 
ter at 1.04 fim at phase 0.40 is 340 Rq, which corresponds to 
29.5 mas (see Fig. 7 in Woodruff et al. "20041. In the present 
calculation, we adopt this value as the diameter of the stel- 
lar disk. Figures]^ andl^) show the predicted visibilities and 
uniform disk diameters, using the HITEMP database and the 
PS97 line list. The visibilities as well as the uniform disk di- 
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Fig. 4. a: Visibilities calculated from the best-fit model for 
o Cet in the 1 1 //m region for a projected baseline length of 
30 m. In all the panels, the green solid lines represent results 
calculated with the HITEMP database, while the blue solid 
lines represent those calculated with the PS97 list, b: Uniform 
disk diameters calculated for a projected baseline length of 
30 m. In a and b, the bandpasses used by WHT03a are marked 
with the dashed lines. The bandpasses used for the observations 
of o Cet at around 896.83 cm"' lie between those shown (see 
WHT03a). The ranges of the uniform disk diameters measured 
in these bandpasses are shown as the hatched regions. The dot- 
ted line represents the photospheric angular diameter adopted 
in the calculation. A flux contribution of 70% from the dust 
shell is included in the calculations of the model visibilities, 
but not in the computation of the uniform disk diameters. Both 
plots (a and b) are redshifted by 0.248 cm ' with respect to 
the rest wavenumber, which corresponds to a radial velocity of 
83 km s"', to match the observed spectrum shown in Fig.|2 c: 
The normalized intensity profile at 896.85 cm ', which is spec- 
trally convolved with the same top-hat function as used in the 
calculation of the spectrally convolved visibility shown in a. d: 
The spectrally convolved squared visibilities at 896.85 cm"', 
plotted together with the visibility points observed by WHT03a 
(filled diamonds), e: The same spectrally convolved visibility 
squared as in d, but at higher spatial frequencies. 



ameters shown in the figure are calculated for a projected base- 
line length of 30 m, which is approximately the mean of the 
baseline lengths used by WOO, WHT03a, and WHT03b. The 
presence of the extended dust shell lowers the visibility by an 
amount equal to the fraction of the flux contribution of the dust 
shell, and therefore, the visibilities resulting from the stellar 
disk and the warm H2O envelope are lowered by a factor of 
0.3 to account for the flux contribution of the dust shell in this 
wavelength region. Note, however, that the uniform disk di- 
ameters are computed from the visibilities excluding the dust 
shell, because the effect of the presence of the dust shell is al- 
ready taken into account in the determination of the uniform 
disk diameters by WOO, WHT03a, and WHT03b. The band- 
passes used in the ISI observations are marked with the dashed 
lines, and the ranges of the uniform disk diameters observed 
within these bandpasses are shown as the hatched regions. 

Figure 0J) demonstrates that the predicted uniform disk di- 
ameters in the bandpasses centered around 896.83 cm"' are in 
agreement with those observed. The ISI observations within 
these bandpasses were carried out at various phases between 
0.92 and 0.36, and therefore, it should be kept in mind that the 
range of the observed angular diameters shown in the figure 
is aff'ected by a variation of the angular diameter with phase. 
The uniform disk diameters predicted within the bandpass cen- 
tered at 895.1 cm"' are systematically lower than the observed 
value of 55.88 + 0.74 mas. However, the ISI observation for 
this bandpass was carried out at phase 0.26, while the high- 
resolution TEXES spectrum of o Cet, which we used for deriv- 
ing the parameters of the water vapor envelope, was obtained 
at phase 0.36. Furthermore, WHT03a note that the uniform 
disk diameter of o Cet measured at 11.1419 //m increased by 
11% (~6 mas) from 2000 to 2001, which they attribute to the 
non-periodic variation in the angular size and/or asymmetries. 
Therefore, the difference between the predicted angular diam- 
eters within the bandpass centered at 895.1 cm"' and that ob- 
served can be attributed to the mismatch in phase as well as this 
non-periodic variation and/or asymmetries of the stellar disk 
and the warm water vapor envelope. 

Figure |^ shows the intensity profiles predicted at 
896.85 cm"', using the HITEMP database (green solid line) 
and the PS97 line list (blue solid line). The spectrally convolved 
squared visibilities at this wavenumber are plotted as a function 
of spatial frequency for both cases with the HITEMP database 
and the PS97 line list in Fig.|3l, together with the observed val- 
ues for o Cet, which were read off Fig. 1 of WHT03a. The fig- 
ure reveals that the predicted squared visibilities are in agree- 
ment with those observed at spatial frequencies higher than 
~ 2.6x 10^ rad"', while the agreement is slightly poorer at lower 
spatial frequencies. However, this small discrepancy can be due 
to the variation of the observed angular diameter mentioned 
above, and therefore, is not regarded as serious disagreement. 
Contemporaneous spectroscopic and interferometric observa- 
tions would be desirable for further constraining our model for 
the warm water vapor envelope. Figure|5 shows the predicted 
squared visibilities at higher spatial frequencies, which corre- 
spond to baselines as long as ~100 m. Mid-infrared interfer- 
ometric observations with such long baselines would also be 
useful for examining our model further 



K. Ohnaka: Warm water vapor envelope in Mira variables 




Fig. 5. Uniform disk diameters around 10.8844 fim (upper 
panel) and 1 1 .0856 fj.m (lower panel) calculated from the best- 
fit model for o Cet for a projected baseline length of 30 m. In 
both panels, the green solid lines represent results calculated 
with the HITEMP database, while the blue solid lines repre- 
sent those calculated with PS97 list. Both plots are redshifted 
by 0.248 cm"' in wavenumber, which corresponds to a radial 
velocity of 83 km s"' , to match the observed spectrum shown in 
Fig. 121 The bandpasses used by WHT03a are marked with the 
dashed lines, and the ranges of the angular diameters measured 
by WHT03a are represented as the hatched regions. The dot- 
ted lines represent the photospheric angular diameter adopted 
in the calculation. 



We perform the same calculation for the two bandpasses 
in the regions around 10.8844 fim and 11.0856 jum. Within 
the bandpasses in these wavelength regions, strong emission 
and absorption features due to H2O are observed, as shown in 
the upper two panels of Fig. |2l Figure [S] shows the uniform 
disk diameters calculated with the HITEMP database and the 
PS97 line list in these wavelength regions. The figure demon- 
strates that the observed uniform disk diameter within the band- 
pass at 918.56 cm"' can be well reproduced by our model, 
while the uniform disk diameters predicted within the band- 
pass at 902.0 cm"' are lower than the observed value. It should 
be noted, however, that the ISI observations for these band- 
passes were carried out at phases diff'erent from the phase 0.36 
at which the TEXES spectrum was obtained: phase 0.95 for 
the bandpass at 918.56 cm"' and phase 0.20 for the bandpass 
at 902.0 cm"'. This mismatch in phase as well as the system- 
atic increase of the angular size from 2000 to 2001 mentioned 
above introduces ambiguities in comparing the predicted and 
observed uniform disk diameters. Given these ambiguities, the 
agreement between the predicted and observed uniform disk di- 
ameters within the two bandpasses containing the strong H2O 
features can also be regarded as fair 
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Fig. 6. Comparison between the ^-band visibilities observed 
for o Cet and those predicted by the best-fit model for this ob- 
ject. The green solid line represents the visibility calculated 
with the HITEMP database, while the blue solid line represents 
that calculated with PS97 list. The observed visibilities include 
the results obtained by Ridgway et al. ( 1992 ) (phase 0.23-0.36, 
red open circles), the speckle data obtained at 0.26 by Woodruff 
et al. (7UU^\ (red filled circles), and the VLTIA^INCI data ob- 
tained at phase 0.26 (red open triangles), 0.40 (purple open 
squares), 1 .40 (black open diamonds), and 1 .47 (orange filled 
diamonds). 

3.1.3. Near-infrared angular diameter 

We have shown that the high-resolution 11 jjm spectrum of 
o Cet and its 1 1 fim angular diameter almost twice as large 
as those measured in the near-infrared can be reproduced by 
our model for the warm water vapor envelope. However, if the 
effect of the warm water vapor envelope on the near-infrared 
angular diameter is as large as that on the 1 1 //m angular diam- 
eter, our model may not reproduce the observed near-infrared 
visibilities and angular diameters of o Cet. If it is the case, 
our model for the warm water vapor envelope as well as the 
photospheric angular diameter of 29.5 mas that we assumed 
based on the ^-band VINCI observations may not be justified. 
Therefore, we calculate the visibility as well as the angular 
diameter of o Cet in the near-infrared using the same best-fit 
model for o Cet and compare with the observed values. 

Although a number of near-infrared interferometric obser- 
vations of o Cet are presented in the literature, only a limited 
amount of data are available up to now for phases near 0.36. 
Given uncertainties in the determination of a zero point of pul- 
sational phase, we can compare with data obtained at phases 



of approximately 0.36 ±0.1. Ridgway et al. (I1992> obtained 
/T-band visibilities between phase 0.23 and 0.36 with base- 
line lengths of 8 and 12 m. The /T-band interferometric data 



of o Cet presented by Woodruff et al. (I2004> include data ob- 
tained at phases 0.26, 0.4, 1 .4, and 1 .47 (phases larger than 1 
mean the next pulsation cycle), together with speckle interfer- 
ometric data obtained with the 6 m telescope at the Special 
Astrophysical Observatory in Russia. The data obtained by 
Ridgway et al. (1992) as well as by Woodruff et al. ( 2004 ) were 
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acquired with a broadband K filter, which means that the de- 
rived apparent size is affected by molecular spectral features 
originating in the warm molecular envelope, in particular, H2O 
lines located at the shorter and longer edges of the K band. We 
calculate the uniform disk diameter in the K band in almost 
the same manner as in the 1 1 yum region. The only difference 
is that we calculate the monochromatic visibility squared from 
the monochromatic intensity profile, and this monochromatic 
visibility squared is spectrally convolved for comparison with 
the above /T-band visibility measurements. We include only 
H2O lines in the calculation, and the ^-band filter is approx- 
imated with a top-hat function centered at 2.2 /im with A/l = 
0.4 jum. The /T-band speckle visibility obtained at phase 0.26 
shows no evidence of an extended dust shell, whose presence 
would result in a steep drop of visibility at low spatial frequen- 
cies (see Fig. 2 in Woodruff et al. .2004j . While a possible ap- 
pearance of significant dust emission in the K band at phases 
other than 0.26 cannot be ruled out, it seems to be unlikely that 
dust emission has a noticeable effect on the /T-band visibility at 
phase 0.36. Therefore, no flux contribution from the dust shell 
is included in the calculation of visibilities in this wavelength 
region. Figure|6lshows a comparison between the observed K- 
band visibilities mentioned above and those predicted by the 
best-fit model for o Cet. The figure illustrates that the observed 
visibilities are well reproduced by our model. 

We also compare the predicted ^-band uniform disk di- 
ameters with that observationally derived. From the predicted 
visibilities shown in Fig. |6l the uniform disk diameter is de- 
rived for a projected baseline length of 15 m, which is roughly 
the mean of the baselines used in the observations presented 
by Woodruff et al. QQQA). The /T-band uniform disk diame- 
ters predicted from the best-fit model for o Cet are 34. 1 mas 
with the HITEMP database and 34.9 mas with the PS97 line 
list. These predicted values are in agreement with the uniform 
disk diameter of 33.27 + 0.33 mas at phase 0.40 obtained by 
Woodruff et al ( 2004 1, and this agreement lends support to our 
model for the warm water vapor envelope as well as the photo- 
spheric angular diameter adopted in the present calculation. 

Ireland et al. ("2004'i measured the angular size of o Cet in 
the near-infrared between 1.2 and 3.6 /^m by aperture mask- 
ing with the Keck telescope. The wavelength dependence of 
the angular size in the near-infrared appears to coincide with 
the wavelength dependence of the opacity of H2O, which sug- 
gests that the observed near-infrared angular size may also be 
explained by our model for the warm water vapor envelope. We 
calculate the angular size, which is the Gaussian FWHM in this 
case, using the best-fit model for o Cet, with only H2O lines 
taken into account. The calculated visibility is convolved with 
a spectral resolution of 120, which approximately corresponds 
to the spectral resolution used for the Keck observations, and 
the Gaussian FWHM at each wavelength is derived for a pro- 
jected baseline length of 10 m. Figure shows a comparison 
between the angular sizes predicted from our best-fit model 
for o Cet from 1.5 to 4.0 yum and that observed by Ireland et 
al. (2004 1. The predicted angular sizes are in good agreement 
with the Keck observation, suggesting that the wavelength de- 
pendence of the angular size in the near-infrared largely re- 
flects the opacity of H2O lines. However, the Keck observation 




Fig. 7. Wavelength dependence of the angular size (Gaussian 
FWHM) predicted from the best-fit model for o Cet for a pro- 
jected baseline length of 10 m. The red dots represent the an- 
gular size of o Cet observed by Ireland et al. (2004). The green 
solid line represents the result obtained with the HITEMP 
database, while the blue solid line represents that obtained with 
the PS97 line list. No flux contribution from the dust shell is in- 
cluded in the calculation shown. 



of o Cet was carried out at phase 0.95, which is far from the 
phase 0.36 of the TEXES observation of o Cet. Therefore, the 
apparent agreement between the observed and predicted near- 
infrared angular sizes of o Cet should be regarded as prelim- 
inary, and this highlights again the importance of coordinated 
spectroscopic and interferometric observations. 



3.2. R Leo 

3.2.1. 11 jum and near-infrared spectra 

For R Leo, we find that the observed 1 1 /im spectrum and uni- 
form disk diameters obtained by WHT03b can be best repro- 
duced by a model with Jhot = 2000 K, 7?hot = 1.7 7?*, A^hot = 
1 X 10^' cm-2, Tcooi = 1200 K, R^ooi = 2.2 R^, and A^cooi = 
7x10^' cm"^. FigureOshows the synthetic spectra in the 1 1 /zm 
region predicted by the best-fit model for R Leo, together with 
the TEXES spectrum obtained at phase 0.75, which was read 
off Fig. 2 in WHT03b. The synthetic spectrum calculated with 
the HITEMP H2O line list is represented with the green solid 
line, while that calculated with the PS97 H2O line list is repre- 
sented with the blue solid line. We estimated the flux contribu- 
tion of the circumstellar dust shell in the 1 1 jjm region based 
on the modeling of the spectral energy distribution (SED) and 
the visibility presented by Schuller et al. (2004 i. Their model 
SED for R Leo, which can well reproduce the photometric ob- 
servations in particular in the infrared, shows that the flux con- 
tribution of the dust shell at 1 1 fim is approximately 70% of the 
total flux (see Fig. 3 of Schuller et al. 2004i, and we adopted 
/dust = 0.7 in our calculation for R Leo. The synthetic spec- 
tra are redshifted by 7 km s"' to match the observation (see 
WHT03b). The figure shows that the synthetic spectra can re- 
produce the observed spectra without prominent spectral fea- 
tures in the bandpasses used for the ISI observations (marked 
with the dashed lines), which is a result of the filling-in effect 
due to the emission of the H2O lines from the extended enve- 
lope. The approximate depths and heights of the observed fine 
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Fig, 8. Spectra of R Leo in the 1 1 fim region at phase 0.75. The 
red dots represent the observed spectrum of R Leo presented 
in WHT03b, while the green and blue solid lines represent the 
calculated spectra using the HITEMP database and the PS97 
line list, respectively. The parameters of the best-fit model for 
R Leo are That = 2000 K, T^hot = L7 R-^, N^^ot = 1 x lO"' cm^^, 
Tcooi = 1200 K, /?cooi = 2.2 R-^, and A^cooi = 7 X lO^^ cm^l The 
synthetic spectra are convolved with a Gaussian with a FWHM 
of 0.013 cm"' to account for the effects of the instrument as 
well as of the macro-turbulent velocity, and are redshifted by 
7 km s ' to match the observation. The dashed lines repre- 
sent the bandpasses used in the ISI observations. The positions 
of the H2O lines whose intensity at 2000 K is stronger than 
1 X 10"^"* cm molecule ' are marked with upper ticks (PS97 
line list) and lower ticks (HITEMP database). These line po- 
sitions are also redshifted by 7 km s ' with respect to the rest 
wavenumber. 



spectral features are also reproduced, as in the case of o Cet 
discussed above. 

We also compare the predicted spectra in the near-infrared 
with that observed. For R Leo, the near-infrared spectrum ob- 
tained at phase 0.5 1 by Strecker et al. J1978> using the Kuiper 
Airborne Observatory (KAO) would be the spectrum closest 
to that expected at phase 0.75, at which the high-resolution 
1 1 /L/m spectrum of R Leo was obtained. Figure|5]shows a com- 
parison between the spectrum of R Leo observed by Strecker 
et al. 61978) and those predicted from the best-fit model for 
R Leo, which are convolved to match the spectral resolution of 
50 used by Strecker et al. ("1978). As we will see below, the K- 
band visibility observed with the Keck telescope (Monnier et 
al. l2004> does not show a steep drop at low spatial frequencies 
characteristic of an object with an extended dust shell. In ad- 
dition, the SED modeling of Schuller et al. ( 2004 ) shows that 
the flux contribution of the dust shell is ~3% at most in the 
near-infrared (A ^ 4 fim). Therefore, the flux contribution of 
the dust shell is neglected in the calculation of the synthetic 
spectra. A glance of Fig. |9] reveals that the observed spectrum 
can be reasonably reproduced by the best-fit model for R Leo 
for phase 0.75, and the difference between the observed and 
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Fig. 9. Spectra of R Leo in the region between 1 .2 and 4.0 fim. 
The red dots represent the KAO spectrum of R Leo observed 
at phase 0.51 by Strecker et al. (1978 1, while the green and 
blue solid lines represent the spectra from the best-fit model 
for phase 0.75, using the HITEMP database and the PS97 Une 
list, respectively. The parameters of the model are given in the 
legend to Fig.|8l The synthetic spectra are convolved to match 
the spectral resolution of the KAO observation (R ^ 50). 




Fig. 10. Uniform disk diameters calculated from the best-fit 
model for R Leo in the 1 1 fj.m region for a projected baseline 
length of 30 m. The green solid line represents the result cal- 
culated with the HITEMP database, while the blue solid line 
represents that calculated with the PS97 line list. The band- 
passes used by WHT03b are marked with the dashed lines. 
The ranges of the diameters measured in these bandpasses are 
shown as the hatched regions. The dotted line represents the 
photospheric angular diameter adopted in the calculation. The 
calculated uniform disk diameters are redshifted by 0.021 cm~' 
with respect to the rest wavenumber, which corresponds to a 
radial velocity of 7 kms ', to match the observed spectrum 
shown in Fig.|S] 

predicted spectra might be attributed to the mismatch in phase 
between the observation and the model. 

3.2.2. Angular diameters measured at 11 //m and in 
tine near-infrared 

We now compare the 1 1 /im angular diameters predicted from 
our model for R Leo with those obtained by the ISI observa- 
tions. Figure[roishows a comparison between the observed uni- 
form disk diameters and those predicted from the best-fit model 
for R Leo. We tentatively adopt a photospheric angular diame- 
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Fig. 11. Comparison between the A'-band visibilities of R Leo 
observed by Monnier et al. ("2004") and those predicted by the 
best-fit model for this object. The red filled circles represent 
the Keck observation, while the red filled diamonds represent 
the IOTA observation. The predicted visibilities which are con- 
volved with a spectral resolution of 43 to match the Keck ob- 
servation are plotted by the green solid line (calculated with 
the HITEMP database) and the blue solid line (calculated with 
the PS97 line list). The dashed lines represent the visibilities 
convolved with a spectral resolution of 7 to match the IOTA 
observation. The result calculated with the HITEMP database 
is plotted by the green dashed line, while that calculated with 
the PS97 line list is plotted by the blue dashed line. 



ter of 28 mas based on the measurements carried out at phase 
0.24-0.28 by Perrin et al. (1999) and examine the validity of 
this assumption by comparing with observed angular sizes in 
the near-infrared below. The uniform disk diameters are cal- 
culated from the visibilities excluding the dust shell for a pro- 
jected baseline length of 30 m, as in the case of o Cet. FigureFTUI 
illustrates that the uniform disk diameters predicted within the 
bandpass at around 896.85 cm"' are in agreement with those 
measured with ISI. However, two measurements within this 
bandpass were carried out at different epochs and phases, and 
therefore, it is necessary to compare the model predictions with 
the observations carried out at similar phases to that at the time 
of the TEXES observation of R Leo. 

The ISI observations were carried out in 1999 October 
(phase 0.43) and in 2001 October and November (phase 0.77) 
for the bandpass at 896.85 cm ' and in 2001 November (phase 
0.81) for the bandpass at 895.1 cm"'. Since our model for 
R Leo is based on the fitting to the TEXES spectrum obtained 
in 2000 December at phase 0.75, the predicted uniform disk 
diameters should be close to the values observed at similar 
phases, that is, 62.62 +1.14 mas at 896.85 cm"' at phase 0.77 
and 64.24 + 1.12 mas at 895.1 cm"' at phase 0.81. Although 
the predicted uniform disk diameters in the two bandpasses 
are somewhat smaller than these observed values obtained at 
phases close to 0.75, this difference may be attributed to a 
cycle-to-cycle variation and/or deviation from circular symme- 
try as WHT03a discuss for o Cet. Given such a possible tempo- 



ral variation of the angular size and a presence of asymmetries, 
we can conclude that our model for the warm water vapor en- 
velope can fairly reproduce the observed 1 1 jum uniform disk 
diameter of R Leo. 

As in the case of o Cet, we compare near-infrared visi- 
bilities predicted by the best-fit model for R Leo with those 
observed at phases close to 0.75. Monnier et al. ( I2004t have 
recently combined interferometric data from Keck aperture 
masking and the Infrared Optical Telescope Array (IOTA) for 
ten objects including R Leo. They observed R Leo at phase 0.7, 
very close to the phase 0.75. Figure^Jshows a comparison be- 
tween the ^-band visibilities obtained by Monnier et al. <2004> 
and those predicted from the best-fit model for R Leo. The ob- 
served visibility at low spatial frequencies (:$ 4 x 10'' rad"') 
does not show a sharp drop, which would be present if there 
were flux contribution from the dust shell. Therefore, we in- 
clude no dust emission in the calculation of the TiT-band visibil- 
ities. The predicted visibilities are spectrally convolved with a 
Gaussian with a FWHM of 0.053 ^^m {R ^ 43) to match the 
width of the narrow-band filter used for the Keck observation. 
For comparison with the IOTA data, which were obtained with 
a broadband K' filter, the predicted visibilities are convolved 
with a top-hat function centered at 2. 16 /im with A/l - 0.32 //m 
(R ^ 7). Figure^Jdemonstrates that the visibility obtained by 
Keck aperture masking (red filled circles) can be well repro- 
duced by the best-fit model, while the predicted visibilities are 
too low compared to the visibility points obtained with IOTA 
(red filled diamonds). This discrepancy at high spatial frequen- 
cies may be attributed to the simplicity of our model and/or a 
cycle-to-cycle variation which might have taken place between 
the Keck+IOTA observations (2000 January and February) and 
the TEXES observation (2000 December). 

Chagnon et al. J2002t measured L'-band visibility of R Leo 
at two epochs, phase ~0.6 and ~0.8, both of which are rather 
close to the phase 0.75. Figure^lshows a comparison between 
the observed L'-band visibilities of R Leo and those predicted 
by the best-fit model for this object. Since the flux contribu- 
tion of the dust shell is expected to be ~3% at most in the L' 
band as mentioned above, no flux contribution from the dust 
shell is included in the calculation of the L'-band visibilities. 
The predicted visibiUties are spectrally convolved with a top- 
hat function centered at 3.8 //m with A/1 = 0.54 //m to match 
the observed data. The figure illustrates that the observed visi- 
bilities can be fairly reproduced by our model, given the slight 
mismatch in phase as well as a possible cycle-to-cycle varia- 
tion of the warm water vapor envelope. This agreement of the 
L'-band visibility lends further support to the photospheric an- 
gular diameter that we adopted and our basic picture of the 
warm water vapor envelope. 

We also compare the wavelength dependence of the angu- 
lar size predicted in the near-infrared with that observed by 
Ireland et al. (2004) . Their observation of R Leo in the wave- 
length region between 1.2 and 3.6 //m with the Keck telescope 
was carried out at phase 0.74, which is very close to the phase 
of 0.75 at which the high-resolution 1 1 yum spectrum was ob- 
tained. We calculate the near-infrared angular size using the 
best-fit model for R Leo in the same manner as in the case of 
o Cet. Comparison between the observed Gaussian FWHM of 
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Fig. 12. Comparison between the L' -band visibilities of R Leo 
observed by Chagnon et al. ("2002^ and those predicted by the 
best-fit model for this object. The black filled circles repre- 
sent the visibilities observed in 2000 March (phase 0.81-0.82), 
while the red filled circles represent those observed in 2000 
November (phase 0.61-0.64). The green solid line represents 
the predicted visibility calculated with the HITEMP database, 
while the blue solid line represents that calculated with the 
PS97 Une list. 
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Fig. 13. Wavelength dependence of the angular size (Gaussian 
FWHM) predicted from the best-fit model for R Leo for a pro- 
jected baseline length of 10 m. The red dots represent the angu- 
lar size of R Leo observed by Ireland et al. (I2004t . The green 
solid line represents the result obtained with the HITEMP 
database, while the blue soUd line represents the result with 
the PS97 line list. No flux contribution from the dust shell is 
included in the calculation shown. 



R Leo and those predicted from our model shown in Fig. E] 
illustrates that our model can reproduce the observed Gaussian 
FWHM from 1.5 to 3.6 yum quite well. This agreement of the 
near-infrared angular size of R Leo suggests again that the H2O 
lines govern the wavelength dependency of the angular size of 
Mira variables in this wavelength regime. 
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Fig. 14. Spectra of ^ Cyg in the 1 1 fim region at phase 0.4. The 
red dots represent the observed spectrum of x Cyg presented 
in WHT03b, while the green and blue solid lines represent the 
calculated spectra using the HITEMP database and the PS97 
line list, respectively. The parameters of the best-fit model for 
X Cyg are Jhot = 2000 K, i^hot = 1.5 R^, A^hot = 1 x lO^' cm-^, 
Tcooi = 1200 K, Rcooi = 2.2 /;*, and A^cooi = 1 x 10^' cm^^. The 
synthetic spectra are convolved with a Gaussian with a FWHM 
of 0.013 cm"' to account for the effects of the instrument as 
well as of the macro-turbulent velocity, and are blueshifted by 
22 km s"' to match the observation. The dashed lines repre- 
sent the bandpasses used in the ISI observations. The positions 
of the H2O lines whose intensity at 2000 K is stronger than 
1 X 10"^"* cm molecule"' are marked with upper ticks (PS97 
line list) and lower ticks (HITEMP database). These line posi- 
tions are also blueshifted by 22 km s"' with respect to the rest 
wavenumber. 

3.3. X Cyg 

3.3.1. 11 /im and near-infrared spectra 

For X Cyg, we find that the observed 1 1 //m spectrum and the 
uniform disk diameter obtained by WHT03b can be best re- 
produced by a model with Thot = 2000 K, /?hot = 1-5 /?*, A'hot 
= 1 X 10^' cm-2, Teooi = 1200 K, /?cooi = 2.2 /?*, and A^cooi 
= 1 X 10^' cm"^. Figure [T^ shows a comparison between the 
TEXES spectrum of ;t' Cyg, which was obtained in 2001 mid- 
June at phase 0.4, and those predicted by the above best-fit 
model for this object. The flux contribution of the circumstel- 
lar dust shell at 11 yum was estimated from the visibility func- 



tion at 1 1 yum obtained by Danchi et al. J 19941 1. which indi- 
cates a steep drop at low spatial frequencies resulting from the 
presence of the extended dust shell. The amount of the visi- 
bility drop, which is equal to the fraction of the flux contri- 
bution of the dust shell, is approximately 0.6, and we adopted 
/dust = 0.6 in our calculation for x Cyg. The synthetic spec- 
tra are blueshifted by 22 km s ' to match the observation (see 
WHT03b). Figure [3]demonstrates that the observed spectrum 
(red dots) can be well reproduced by the model. Like in the 
cases of o Cet and R Leo discussed above, the fiUing-in effect 
due to H2O line emission from the extended envelope leads to 
the featureless, continuum-like spectra. 

The presence of such a dense water vapor envelope around 
X Cyg, though the H2O column density is smaller than that 
in o Cet and R Leo, may appear to be in conflict with the 
classification of this object as an S star S stars have C/O ra- 
tios close to unity, and therefore, neither oxygen-bearing nor 
carbon-bearing molecules are abundant in the atmosphere, with 
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most of carbon and oxygen atoms locked up in CO molecules. 
In fact, the H2O lines observed for x Cyg in the 4 jjm re- 
gion are significantly weaker than those observed in oxygen- 
rich Miras (Lebzelter et al. 1200 1> . However, there are some 
pieces of observational evidence for the presence of a signif- 
icant amount of water vapor in x Cyg near minimum light. 
Based on the H2O absorption features in the 9000 A region, 
Spinrad & Vardya J1966t estimated the H2O column density 
in X Cyg to be 3 x 10^^ cm"^ at minimum light, which is in 
agreement with the result we obtained above. Wallace & Hinkle 
<1997> present a series of ^-band spectra of x Cyg observed 
at various phases with a spectral resolution of ~3000 using a 
Fourier Transform Spectrometer (FTS), and the spectra exhibit 
an increasing absorption in the region between ~4700 cm ' 
and 4950 cm"' from maximum light (phase 0.00) to minimum 
light (phase 0.42-0.48). Figure [Tst shows the spectra acquired 
at phase 0.00 and 0.42', and the thick red solid line in Fig. 115b 
represents the ratio of these two spectra (i.e., spectrum at 0.42 
divided by spectrum at 0.00). This ratioed spectrum clearly 
shows the enhanced absorption which appears between ~4700 
and 4950 cm ' near minimum light. The thin green solid line 
in Fig. 115b represents the H2O spectrum predicted by the above 
best-fit model for;if Cyg using the HITEMP database. Dyck et 
al. 11984 ) estimated that more than 90% of the observed flux 
between 2.2 and 4.8 jum comes from the star based on the mod- 
eling of Rowan-Robinson & Harris (i983i|. Therefore, the flux 
contribution of the dust shell is not included in the calculation 
of the synthetic spectrum. As the figure shows, the absorption 
observed between ~4700 and 4950 cm ' is reproduced by the 
synthetic spectrum, and the individual absorption features seen 
in the ratioed spectrum can be identified as H2O lines, whose 
line positions are marked with the ticks in the figure. 

Finally, more definitive evidence for the presence of a large 
amount of water vapor in x Cyg near minimum light can be 
found in the near-infrared spectrum obtained with ISO/SWS. 
Vandenbussche et al. (2002 1 present the ISO/SWS spectrum of 
X Cyg obtained on 1998 May 10 (UT 07:39:48) at phase 0.5 in 
the wavelength region between 2.3 and 4.1 yum. We retrieved 
this ISO/SWS spectrum of x Cyg from the ISO data archive 
and plot it in Fig.[^ The spectrum clearly exhibits the absorp- 
tion feature centered at 2.7 jjm due to the H2O vi and V3 fun- 
damental bands. Also plotted in the figure are the spectra pre- 
dicted from the above best-fit model for x Cyg, and the figure 
demonstrates that the predicted spectra are in fair agreement 
with that observed, given the simplicity of our model as well 
as a slight mismatch in phase between the TEXES observation 
and the ISO/SWS observation. 

The above observational results show that water vapor is 
present in the S-type Mira variable ;t' Cyg near minimum light, 
and that the amount of water vapor is comparable to those 
found in the oxygen-rich Miras o Cet and R Leo. In general, 
Mira variables exhibit the strongest molecular absorption near 
minimum light (e.g., Lan9on & Wood 2000), which implies 
that the temperature may be the lowest at minimum light and 
that the formation of water vapor may be significantly pro- 
moted at such low temperatures. However, detailed calculations 
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Fig. 15. Spectra of x Cyg in the region between 4600 and 
4950 cm"', a: The dotted line represents the spectrum of;^^ Cyg 
obtained at phase 0.00, while the solid line represents the spec- 
trum at phase 0.42 (Wallace & Hinkle , 1997 ). b: The thick red 
solid line represents the ratio between the spectra observed at 
phase 0.00 and 0.42. The thin green solid line represents the 
spectrum predicted by the best-fit model for;^' Cyg, whose pa- 
rameters are given in the legend to Fig.^^ The synthetic spec- 
trum is convolved to match the spectral resolution of R ^ 3000 
used in the observations by Wallace & Hinkle tl997j . 



using dynamical atmospheres are necessary for understanding 
the formation of H2O in Mira variables with C/O ratios very 
close to unity. 

3.3.2. Angular diameters measured at 11 /^m and in 
the near-infrared 

In Fig. [21 we compare the angular diameter of;^' Cyg measured 
in the 1 1 /im region and those predicted by the above best-fit 
model for this object. The phase of ;i'Cyg at the time of the 
ISI observation is 0.51 (2001 July and August, see WHT03b), 
which is quite close to the phase 0.4 at which the TEXES spec- 
trum was obtained. We adopt a photospheric angular diame- 
ter of 23 mas based on a /T'-band uniform disk diameter mea- 



Available at_ftp://ftp.noao.edu/catalogs/medresIR/K_band/_ 



sured by Mennesson et al. (I2002> at phase 0.38, which is also 
close to the phase 0.4. As in the cases of o Cet and R Leo, we 
check the validity of the assumed photospheric angular diam- 
eter by calculating the near-infrared angular diameters below. 
The uniform disk diameters are calculated from the visibilities 
excluding the dust shell for a projected baseline length of 30 m. 
Figure [n] shows that the uniform disk diameters predicted in 
the ISI bandpass are in agreement with the observed value of 
39.38 + 4.02 mas. 
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Fig. 16. Spectra of ;^' Cyg in the region between 2.3 and 4.0 //m. 
The red solid line represents the spectrum of ;t' Cyg obtained at 
phase 0.5 with ISO/SWS (Vandenbussche et al. l2002t . while 
the green and blue solid lines represent the spectra from the 
best-fit model for phase 0.4, using the HITEMP database and 
the PS97 line list, respectively. The parameters of the model 
are given in the legend to Fig. ^^ The synthetic spectra are 
convolved to match the spectral resolution of the ISO/SWS ob- 
servation {R ^ 1500), and no flux contribution from the dust 
shell is included. 




Fig. 17. Uniform disk diameters calculated from the best-fit 
model for x Cyg in the 1 1 yum region for a projected baseline 
length of 30 m. The green solid line represents the result ob- 
tained with the HITEMP database, while the blue solid line 
represents the result with the PS97 line list. The bandpass used 
by WHT03b is marked with the dashed lines. The range of the 
diameters measured in this bandpass is shown as the hatched 
region. The dotted line represents the photospheric angular di- 
ameter adopted in the calculation. The calculated uniform disk 
diameters are blueshifted by 0.0658 cm"' with respect to the 
rest wavenumber, which corresponds to a radial velocity of 
-22 km s ', to match the observed spectrum shown in Fig. 1141 



We also compare the angular diameters predicted in the 
near-infrared with those observed. Mennesson et al. (2002 1 ob- 
served X Cyg in the K' and U bands and derived uniform disk 
diameters of 23.24 + 0.08 mas and 30.40;^^2g mas, respectively. 
We calculate the uniform disk diameters in these bands, with 
the K'- and L'-band filters represented with top-hat functions 
centered at 2. 16 yum with AA = 0.32 yum and at 3.8 fim with AA 
= 0.54 jum, respectively. The ^'-band uniform disk diameters 
calculated with the HITEMP database and the PS97 line list for 
a projected baseline length of 20 m are 24.4 mas and 24.7 mas, 
respectively. These predicted /T'-band uniform disk diameters 
are in reasonable agreement with the observed value, which 
was measured at phase 0.38, very close to the phase 0.4 of the 
TEXES observation. The L'-band uniform disk diameters cal- 
culated with the HITEMP database and the PS97 line list ai-e 
27.1 mas and 28.1 mas, respectively. While these predicted L'- 
band uniform disk diameters are also in agreement with the 
observed one within the error of the measurement, it should be 
noted that the L'-band observations of Mennesson et al. (2002) 
were carried out at phase 0.81, which makes it difficult to com- 
pare the observed angular diameter with that predicted by the 
best-fit model based on the TEXES spectrum obtained at phase 
0.4. Therefore, we can conclude that the near-infrared angu- 
lar diameters predicted by our model for x Cyg are in rough 
agreement with those observed, while more rigorous tests are 
necessary using contemporaneous interferometric observations 
in the near-infrared and in the mid-infrared. 



4. Discussion 

Our two-layer model for the warm water vapor envelope around 
three Mira variables can reproduce the observed increase of 
the angular diameter from the near-infrared to the 1 1 yum re- 
gion and, simultaneously, the observed high-resolution 1 1 yum 



spectra which appear to be featureless. It is obvious, however, 
that our ad hoc model is not a unique solution, and that physi- 
cal processes responsible for the formation of such warm water 
vapor layers also remain to be theoretically understood. 

Intensity profiles predicted by dynamical model atmo- 
spheres for Mira variables exhibit structures consisting of two 
components. Tej et al. J2003t show that the formation of H2O 
layers behind shock fronts leads to intensity profiles consist- 
ing of a step-like structure close to the star and/or a tenuous 
halo-Uke structure extending to a few stellar radii. The inten- 
sity profiles predicted by our models for the three Mira vari- 
ables studied here also show such a step-like structure and a 
more extended component, as can be seen in Fig.|4l;. This im- 
plies the possibility that the large-amplitude pulsation in Mira 
variables is responsible for the warm molecular envelope in 
these objects, and the next logical step would be to calculate 
synthetic spectra as well as intensity profiles using dynamical 
model atmospheres such as presented in Tej et al. (2003 it and 
extensive H2O Une lists. 

Hinkle & Barnes (19791 analyzed H2O lines in high- 
resolution near-infrared spectra (4000-6700 cm ') of R Leo 
and identified two distinct H2O layers, that is, a warm com- 
ponent with a temperature of ~1700 K and a cool compo- 
nent with a temperature of ~ 1100 K. This result seems to be 
in agreement with the temperatures of the H2O layers we ob- 
tained for R Leo in the present work. The H2O column densi- 
ties of the cool component derived by Hinkle & Barnes ( I1979> 
range from lO^'-lO^^ cm"^, and these values also agree with 
the 7 X 10^' cm"^ we derived for the cool H2O layer. On the 
other hand, the H2O column densities of the warm compo- 
nent they derived are lO'^-lO^'^ cm"^, which are significantly 
smaller than the 1 x 10^' cm"^ derived here. Hinkle & Barnes 
(1979) also note that the H2O lines originating in the warm 
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component are never predominant compared to those originat- 
ing in the cool component, which shows a marked contrast to 
our result that the hot H2O layer contributes significantly in the 
spectrum as well as in the intensity profile. However, Hinkle 
& Barnes ( I1979> derived the H2O column densities at phases 
0.94, 0.00, 0. 14, and 0.20, but not at phases between 0.5 and 0.9 
because of the blend of H2O lines originating in the two com- 
ponents. Since our model for R Leo is based on the TEXES 
spectrum obtained at phase 0.75, the above discrepancy of the 
H2O column density of the warm component might be due to 
the mismatch in phase. Comparison with synthetic spectra us- 
ing dynamical model atmospheres would be necessary to derive 
physical parameters such as excitation temperature and column 
density from the high-resolution spectra observed near mini- 
mum light which are severely affected by the blend of lines. 



In order to test our model for the warm water vapor enve- 
lope further, contemporaneous spectroscopic and interferomet- 
ric observations are indispensable, given the variations of the 
spectra and the angular size not only with phase but also on a 
time scale longer than the variability period. ISl observations 
using many more bandpasses in the 1 1 fim region and/or with 
even higher spectral resolution would also be useful for this 
purpose. 
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5. Concluding remarks 

We have shown that our two-layer model for the warm wa- 
ter vapor envelope around the Mira variables o Cet, R Leo, 
and X Cyg can reproduce the observed increase of the angu- 
lar diameter from the near-infrared to the 1 1 fim region and, 
simultaneously, the observed high-resolution 11 /im spectra. 
While a number of H2O pure rotation lines are present in the 
wavelength region observed with ISI and strong absorption can 
be expected, the absorption lines are filled in by the emission 
lines originating in the extended part of the water vapor en- 
velope, which leads to the featureless, continuum-like spectra. 
This filling-in effect masks the spectroscopic fingerprints of the 
warm water vapor envelope, but its presence manifests itself 
as an increase of the angular diameter from the near-infrared 
to the mid-infrared: invisible for spectrometers but not for in- 
terferometers. The radii, temperatures, and H2O column den- 
sities of the hot H2O layer in the three Mira variables stud- 
ied here are derived to be 1.5-1.7 R^, 1800-2000 K, and (1- 
5) X 10^' cm"^, respectively. The cool H2O layer is found to 
have temperatures of 1200-1400 K, extending to 2.2-2.5 7?* 
with H2O column densities of (1-7) x 10^' cm"^. Our models 
which reproduce the spectra and angular diameters observed at 
1 1 fim have turned out to also reproduce the spectra and the vis- 
ibilities as well as the angular diameters observed in the near- 
infrared. Comparison between the near-infrared angular sizes 
predicted for o Cet as well as for R Leo and those observed 
with the Keck telescope suggests that the wavelength depen- 
dence of the angular size of Mira variables in the near-infrared 
largely reflects the H2O opacity. 

We have also found evidence for the presence of a large 
amount of water vapor in;(f Cyg near minimum light in the FTS 
spectra obtained in the K band as well as in the ISO/S WS spec- 
trum between 2.3 and 4.1 yum. The observed H2O absorption 
features can be reasonably well reproduced by our water vapor 
envelope model whose parameters are derived from the com- 
parison of the spectrum and the angular diameter observed at 
1 1 yum. While S-type Miras are known to exhibit significantly 
weak spectral features due to oxygen-bearing molecules com- 
pared to oxygen-rich Miras, water vapor of an amount com- 
parable to that in oxygen-rich Miras is present in x Cyg near 
minimum light. 
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